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I. INTRODUCTION 



Supersymmetric quantum mechanics(SUSY QM) has provided a possibility to solve 
analytically some non-relativistic quantum systems. The simplest model in this framework 
was proposed by Witten[l] in the early eights of the last Century After this pioneered work, 
the methods of SUSY QM have quickly developed and some applications arised. Some of 



these applications include the problems related with the superpotential 



, yj , whose gener- 



alization to the case of a matrix superpotential was done a longtime ago. The SUSY QM 
has also applications in the non-relativistic quantum context (jj-p]] and in the case involving 
two or three fields in (l+l)-dimensional model. Others applications include some results 

n n 

connected with self-adjoint extensions [8( and superpotential matrix [9], among others. 



(3 



The classical co nfig urations with domain wall solutions are bidimensional structures in 



l)-dimensions [10|-jl4|. 
15( and Prasad- Sommerfield 



'hey are static, non-singular, classically stable Bogomol'nyi 



161 ] (BPS) soliton (defect) configurations, with finite localized 



energy associated with a real scalar field potential model. 

The BPS states are classical configurations that satisfy first and second order differential 
equations. In a context that stresses the connection with BPS-bound states^, Q, domain 
walls have been exploited. 

Recently, the stability and metamorphosis of BPS states have been investigated [1 
using the framework of SUSY QM, with a detailed analysis of a 2-dimensional N = 
2— Wess-Zumino model in terms of two chiral superfields and composite dyons in N = 2 
supersymmetric gauge theories 20]. Also, the superpotential associated with the linear clas 



sical stability of the static solutions 
were discussed in the literature 



or systems with one real scalar field in (l+l)-dimensions 
3J. However, for solitons associated with three coupled 
scalar fields there are no general rules for finding analytic solutions since the nonlinearity in 
the potential leads to an increasing of the difficulties to solve the BPS and field equations. 

This paper is organized as follows: In Section II, we discuss SUSY algebra with topological 
charge. In Section III, we consider a SUSY model for two coupled scalar fields. In Section 
IV, we present the BPS configurations for three coupled scalar fields. In Section V, we 
define the BPS mass bound of the energy and discuss the stability of BPS states. The 
Schrodinger-like equation and also the Hessian matrix are obtained. In Section VI, a matrix 
superpotential with three-component wave functions is obtained. In Section VII, a specific 



2 



potential model is investigated. Our conclusions are presented in Section VIII. 

II. SUPERSYMMETRY ALGEBRA WITH TOPOLOGICAL CHARGE 

Let us start with a discussion concerning central charges due to the fact that these 
quantities characterize SUSY. To do this, consider the potential model with one single real 
scalar field <j), whose Lagragian is given by 



where U(4>) a well defined continuous function and the Majorana spinor, ip(z), is given by 



A = J d 2 z l - {{d^f + [^y v - u\4>) - WWW] } 



(i) 




(2) 



In this case, the conserved SUSY current can be written as 



(3) 



Therefore, the topological SUSY charge is given by 




(4) 



and, then, we can write 




(5) 




(6) 



In (l+l)-dimensions the SUSY algebra becomes 



Ql = P + = P + P u Q\ = P_ = P -P 1 



(7) 



and 







(8) 
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where 

[i>-(y)^-(x)}+ = S(y-x), [ip + (y),ip + (x)] + = S(y-x), 



[il>-(x),il> + (x)] + = 0. 

In a field theory without soliton solutions these equations are satisfied. However, in a 
field theory with soliton solutions these equations are not satisfied because the surface terms 
for a soliton solution are different from zero, and as a consequence 

r+oo Q 

Q+Q- + Q-Q + = y_ oo dz—{2T{4>)), (9) 

with the superpotential satisfying the relation 

r W -±r-u W . (io) 

Note that the right hand side of Eq. (jHJ) is a scalar, which corresponds exactly to the central 
charge. Thus, the Bogomol'nyi classical bound, for a single particle with mass m , at rest, 
which means that, P + = P_ = m , becomes 

W-rMi, (ii) 

where Mi and Mj represent the vacuum states. It is worth calling attention to the fact that 
this inequality remains valid for soliton and antisoliton solutions at one-loop order. 

III. SUSY FROM TWO COUPLED SCALAR FIELDS 

Let us write the potential V{4>j) i n t ne following SUSY form, analogous to the case 
with one single field only, 

+ Uii^^Uii^fc). (12) 

Thus, the N = 1 algebra can be discussed by investigating the SUSY Lagrangian density in 
(l+l)-dimensions with the following form 

£ = I {d,4>if + I (d^ 2 ) 2 + i& (irfdr - ^ ^ 

4 



m > 



f + OQ Q 

J-oo OZ 



-\w^ 2 ' Iw^ 1 (13) 

where ipi and ifj 2 are Majorama spinors. In this framework, the SUSY current is given by 

s» = (dM-f-ffa + iu^j)^ + {dp<h)i p -fth + ^2(^)7^2, (14) 

and therefore, the conserved supercharges can be expressed as 

Q± = J dz {(do4>i ± di4>i)ip± =F Ui((f)j)ip^} 

+-L J dz {(«9 O 2 ± =F U 2 {ff>j)^} . (15) 

On the other hand, the superpotential W(<f>j) satisfy 

Wi = 1 wr w ' ] (16) 

which leads to the value for a Bogomol'nyi minimum energy 

IV. CONFIGURATIOS WITH THREE COUPLED SCALAR FIELDS 

In this section, we consider classical soliton solutions with three coupled real scalar 
fields, 4>j, (j — 1, 2, 3), in (l+l)-dimensions included in bosonic sector and explain the equal- 
ity of topological and central charges, ipi = 0. The soliton solutions are static, nonsingular, 
classically stable and finite localized energy solutions of the field equations. The Lagrangian 
density for such nonlinear system in the natural system of units (c = h — 1), in a (1+1)- 
dimensional space-time, with Lorentz invariance, is written as 

£(<l>i,dM = li:(dM 2 -v(<i> j ), (17) 

Z 3=1 

where <9 M = z^ = (t,z) with \x = 0,1, <j>j = (pj(t,z) and rf v = diag(+,—) is the 

metric tensor. Here, the potential V(<f>j) = V((j) 1 , (f> 2 , 03 ) is a positive definite function of 4>j. 
The general classical configurations obey the following equation 
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d 2 d 2 d 

W*>-e*+> + ej; v = ' (18) 

which, for static soliton solutions, is equivalent to the following system of nonlinear second 
order differential equations 







V, 07 = 1,2,3), (19) 

where prime denotes differentiation with respect to the space variable. 

There is in literature, a trial orbit method for finding static solutions of Eq. (jT9l) . for 
certain positive potentials. This constitutes what is termed the "trial and error" technique 



Uj. This method has many limitations, notably the need to choose trial orbits. Solutions 
had to be obtained by ingenuity combination with trial and error rather than by systematic 
derivation from the field equations. In this paper we will use the trial orbit method for the 
first order differential equations associated with three real scalar fields, differently from what 



was done by Rajaraman 



111 ], who applied this method to the equation of motion. 



Let us assume that the trial orbit is given by 



G(fa,fa,fa) = 0. (20) 

Thus, we have 

A 3 Q/n 

£g(0i,&,&) = £^ = °- (21) 

Taking Eqs. (TSUj) and (|21|) into account, we can get the constant coefficients in such trial 
orbit by substitution of the vacuum and the BPS states into these equations. 

V. LINEAR STABILITY 

Since the potential V((j)j) is positive, it can be written in the square form analogous to 
the case in which we have just one sing le field [a], as 

1 3 8W 

V(fa) = V(fa, fa, fa) = -J2 Uftti, fa> fa), Ujfa, 2 , fa) = — -, (22) 

i j =1 ocpj 

where W is the superpotential associated with the three fields. 
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Therefore, we can write the total energy given by 



E 



-L 



+oc I 

dz- 

OO A 



■A) 2 + 



+ ( ( p' 3 ) 2 + 2V( ( p,x) 



(23) 



in the BPS form, which consists of a sum of squares and surface terms, as 



E 



(24) 



2' " ' " 2' ! " ' dz 
Note that the first three terms are always positive and thus, the lower bound of the energy 
is given by the fourth term, which means that 



E > 



d 

dz—W[ ( f> 1 (z),4>2(z),Uz)} 



(25) 



where the superpotential W = W[(j)i(z), $2(2), </>3 (z)] will be discussed in what follows. The 
BPS mass bound of the energy which results in a topological charge is given by 



E BPS = T ij = \W[M j ]-W[M i ]\, 



(26) 



where and Mj represent the BPS vacuum states and are the extrema of W. In this case 
the BPS states satisfy the following set of first order differential equations 

0;. = ^(0 1)( /) 2 ,03). (27) 

Now, let us analyze the classical stability of the soliton solutions in this nonlinear system, 
taking into consideration small perturbations around (f)j(z)(j = 1,2,3), namely, rjj. Thus, 
we can write the classical solution of the system as 

fait, z) = fa{z) + Vj (t, z), (j = 1, 2, 3). (28) 

We can expand the fluctuations r)j(t,z) in terms of the normal modes, in the following 

way 



(29) 



where e^ n are real constant coefficients. Thus, the stability equation for the fields turns into 
a Schrodinger-like equation for a three-component eigenfunction \l/ n , 
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where 



n^ n = u 2 J> n , n= 0,1,2,-.., 



(30) 



dz 2 ^ d$[ V 



H 



-V 



-V \ 



a 2 



\ 92 V 



~dz 2 



v F (z) 



(31) 



a 2 y -JL+^v I 

d<t>zdcj>2 dz 2 d(K / 



\<Pj=<t>j( z ) 



with I being the (3x3)-dentity matrix and Vp(z) the (3x3) fluctuation Hessian matrix. The 
excited modes are, thus, given by 



1 !&,„(*) ^ 



(32) 



V vsA z ) ) 

Since Vp(z) is a symmetric matrix and Ti is Hermitian, thus the eigenvalues u\ of 7i are 
real. 

The Schrodinger-like equation (|30|) and the Hessian matrix Vp(z) in Eq. (I3TI) are obtained 
by taking a Taylor expansion of the potential V((j)j) in terms of r]j(t,z) and retaining the 
first order terms in the equations of motion. 



VI. POTENTIAL MODEL WITH THREE SCALAR FIELDS 



As an application of this formalism, let us consider the following potential 



1 / TV? 

V= V(<f> 1 ,(j) 2 ,(i> 3 ) = -{\tj)l + a ( f>l + a<f>l- — 



1 2 

+ - (-«0i0 2 + Ml - fo) 



(33) 



where a > and (3i > 0. This is a generalized potential for three scalar fields which was 
constructed from the potential discussed recently 18], for two scalar fields. Note that the 
symmetry Z 2 xZ 2 is preserved only if </> 2 — or if (3\ = (3 2 — 0. 



S 



The corresponding superpotential in a field theory model is given by 



W^j) = - ^<j>\ - aMl ' f Ml + (hMl ~ Pzfo + (34) 

It is required that 4>j, satisfy the BPS state conditions, which are expressed by the 
following equations 



o 

Tfi 

*i = -^4>\ - «02 + — - a 4>l 

b' 2 = -2a<p l( j) 2 + (3 2 4>l - p 2 



1 - 3 (-«0i + 2/? 2 2 + aft) (35) 



and the superpotential W^(0j) satisfy |^ = C/j(j = 1,2,3). 

Note that the BPS states saturate the lower bound, so that Ebps — \Wij\ is the cen- 
tral charge of the realization of N — 1 SUSY in (l+l)-dimensions. Thus, the vacua are 
determined by the extrema of the superpotential. Therefore, the condition 

dW 

Wj =«, 3 — 1-2,3 (36) 

provides the vacuum states Mj = (</>rv, 4>2v, 03v) whose values must satisfy the following 
equations 

Tf? 1 

-\(t>\ - a(j) 2 2 + — - -a03 = 
-2a0i02 + - ft = 

+ 2ft0 2 + aft = 0. (37) 

In order to obtain an explicit form of the vacuum states, let us consider the cases (p 2 v = 
03v = and </> lv = 4>3v = 0, respectively. Thus, we obtain the four vacuum states, which 
are given by 



(m \ 
M 2 =(y,0,0j 
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.U, = I 0. -///\ — .0 




M 4 



0,7m/ — ,0 
Aa 



(38) 



It is easy to verify that these vacuum states are satisfied by the equations given in (1371) . 
for /?2 = and aX > 0. Now, let consider 0rv = Pi and 02v = in Eq.( |37j) . In this case, we 
obtain two additional vacuum states, which are 



Ms 



M fi 



Pi A 



\ 



-{—- m 

a \ A , 



(39) 



for /? 2 = 0, — y < /3i < y an d a ^ 0. Therefore, in this case the components of the tension 



are 



12 



4 m 3 
T21 = 3 



Tl3 


= r 31 


— ^24 


— T42 


= T 41 


= T 14 


— ^23 — 


Tl5 


= r B1 


= T 25 


= T 52 


= T 16 


= T 61 


— T 2 q = 


T34 


= ?43 


= = 


^56 = 


^65 




- 


T 35 


= T 53 


= T45 


= ^54 


= 2~36 


— ^63 



T32 — -^r 
Tq2 



2 m 3 

3X 2 
2 m 3 



A 

2' 



m 
T 



1 • 



(40) 



From the results given by Eq. fj40l . we see that the potential presents two non-topological 
sectors, which are non-BPS sectors, namely, T 34 and T 56 , and twelve BPS topological sectors. 
Now, let us specialize to the (3x3)-matrix superpotential, W, with (3 2 = 0, which is given 

by 



W 



\ 



(41) 



2A0i a<p2 cfp'i 
acf)2 ctcf)i 

This superpotential satisfies the Ricatti equation associated with the non-diagonal fluctua- 
tion Hessian matrix, Vf(z), which is written as 
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W 2 + W = V F (z) 



V Fn (z) V F12 (z) V F13 (z) 



\ 



V F12 (z) V F22 {z) V F23 {z) 
\V F13 (z) V F23 (z) V F33 (z) J 



(42) 



I0=<K 2 )>X=xO) 

where the elements of V F (z), denoted by Vpij(z),aie given by the following relations 



V FU = 6A 2 2 + a 2 (40 2 + 2 ) + 2A (a<& + l -a<g - ^- 



V F22 = Qa 2 <pj + 2a(2a + A)0 2 + (a 2 + 4/3 2 )^ 2 - 2a— 



m 

V F12 = V 21 = 4Aa0i0 2 + 8a 2 0i0 2 - 2a(2/? 2 2 - (3 2 ) 



Vfi3 = V31 = 2a(A + a)0i03 - 4a/3 2 02^3 - 2a 2 /? 



1<?3 



V F33 = -a 2 4>\ + 6/3|0| - 2&(2a0i0a + fa) + (-o;0i + 2/? 2 02 + a^) 2 
^32 = V F23 = 2a 2 2 03 + 40 3 /5 2 (-a0i + 2&0 2 + "A)- 



(43) 



The Ricatti equation (j4"2l only holds for the BPS states. According to Witten's modeljlj, [4j , 
we have 



A ± = ±1— + W(z 



(n) 

SUSY 



» 



(44) 



^ H (^) 

where ty^\z) are three-component eigenfunctions. In this case, the graded Lie algebra of 
the SUSY QM for the BPS states may be realized as 



1x6 



SUSY 



[Q-,Q- 



A + A- 
A-A + 



6X6 



H- 
H 4 



[Hsusy,Q±L = = (Q^) 2 = (Q + ) 2 , 
where Q± are the (6x6) supercharges of the Witten model and is given by 



(45) 



(46) 



q_ = _ ® A , Q+ = Ql 



A+ 




0+ ®^ 



(47) 



with the intertwining operators, A ± , given in terms of (3x3)-matrix superpotential, Eq. fj44l) . 
and a± = ~(cti ± icr 2 ), with o\ and a 2 being Pauli matrices. Note that the bosonic sector of 
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Hsusy is exactly the fluctuating operator given by H- = H = -I^r + Vf{z), where V = 
\~f(z) is the non-diagonal fluctuation Hessian matrix. The supersymmetric fluctuation 
partner operator of H- is 



H+ = A~A + = A + A- + [A~, A + ] =H-- W'(z), 
so that the SUSY partner is given by V + = V_ — ~W'(z). 



(48) 



The Ricatti equation given by (1421) is reduced to a set of first-order coupled differential 
equations. In this case, the superpotential is not necessarily defined as W(z) = ~^to)~^' l P^-\ z )i 
as in the case of a system described by a one-component wave function in the framework of 



SUSY QM 



Therefore, as the zero-mode is associated with a three-component eigenfunction, : i 
one may 



^¥°\z)=w¥°\z) 

dz 



(49) 



from which we find the following zero mode eigenfunction 



(0) 



\ 



UMi) , (50) 

V W j 

where Ui (z = 1,2,3) are given by the BPS states (1351) . 

Now, let us show that the uj n 21 s are non- negative. To do this, consider the bilinear form 
of 7i given by 



U = A + A~ 



(51) 



where 



A' 



(A 



a 2 



(52) 



a x Ji 12 Jii 3 \ <> ■ i 

A 2 \ &2 "4-23 = 
V ^31 ^32 «3 / V a J ) 

with the obvious identification of the elements of TZ(<pi) and the following expressions for 
;he operators that appear in the analysis of classical stability associated with a single field 
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a 3 



~+2\<j> 1 , 
dz 

d 

-— + afa, 
dz 

—— + afa — afc, (53) 



where 

Ai 2 = afa = A21, A n = a<f) 3 = A 31 , ^23 = = ^32- 
Since a+ = (ajY and hence A + = (Ay, we find 

W)« = -& + & (54) 

which are exactly the diagonal elements of 7i. It is worth calling attention to the fact that 
the linear stability is satisfied, which means that 

uj 2 n =< H >=< A + A- >= (A~^ n )\A~^ n ) > 0, (55) 
and therefore the u n 2y s are non-negative. 

VII. PROJECTIONS ON SCALAR FIELDS 

Let us now consider a projection on the (fa, fa) plane in order to find an explicit form 
of domain walls using the trial orbit method. In this case, if we choose fa = in Eq. ( 1331) 
and the following trial orbit 

G(0i,0 2 ) = Cl ^ + c 2 ^ + C3 = o, (56) 

we get from g = j^fa + j^fa = and using the BPS states (j33L that c x = 1, c 2 = 

2 



2 

and C3 = Thus, the resulting elliptical orbit is 



a , m 2 



or 



* + (i=ao* = 7?- < 57 > 



13 



for a = j. These provide the following BPS solutions 



^3 = 0, 



m . tyi . 
— tanh(—z) 

±V2^sech&) 
A z 



(59) 



which connect the vacua (y,0,0) and (— y,0,0). Note that 



i m 

?1 — > ±— , for z — > ±oo 
A 



0, for z — > ±oo. 



This result corresponds to the same one obtained recently 18 
two scalar fields are taken into account. 



for BPS solutions when 



In this situation, the matrix superpotential in SUSY QM, W(z), becomes 



V 



4tanh(fz) ±V2sec/i(f^) 
± v / 2sech(f2) tanh(fz) 
4 



(60) 



where we have used the BPS states in terms of a projection on the (0, x) plane. This provides 
the following bosonic zero-mode 



d 

dz 



,(o) 



(0), 



( sech 2 (fz) 
±v / 2tanh(f z)sech(f z) 




V 



for a — j, where N is the normalization constant. However, in Ref. 
A and a is given by a — 



/ 



(61) 



181 ] the relation between 



VIII. CONCLUSIONS 



In this paper, we considered the classical stability analysis for BPS domain walls associ- 
ated with a potential model of three coupled real scalar fields, which obeys the non-ordinary 
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super symmetry (SUSY). The approach of effective quantum mechanics provides a realization 
of SUSY algebra in the three-domain wall sector of the non-relativistic formalism. 

The components of the tension given in f!4Up were deduced from the charge central prop- 
erties in the model that present N — 1 SUSY. From a three-field specific potential model we 
found two null tensions which correspond to non-topological sectors, and other topological 
sectors, which depend on the manifold of vacuum states, Ty = \W[Mj] — W[Mi}\, where Mi 
and Mj represent the vacuum states. 

We have shown that the positive potentials with a square form lead to three-component 
non-negative normal modes u n 2 > 0, analogous to the case with a single field |2J, so that 
the linear stability of the Schrodinger-like equations is ensured. 

We have seen that domain walls associated with the three-field potentials have features 
that are not present in the one-field models. The BPS states which connect the vacua exist 
within the stability domain and minimize the energy. Thus, they provide a realization of 
the supersymmetric quantum mechanical algebra for three-component eigenf unctions. From 
the stability equation, we have found an expression for the matrix superpotential, satisfying 
the Ricatti equation, within the context of unidimensional quantum mechanics. 

We also deduced an (3x3)-matrix explicit form of the SUSY QM superpotential from a 
field-theoretic superpotential model in (1 + 1) -dimensions. A general three-component zero- 
mode eigenfunction is deduced, but its explicit form is found only for the projection on the 
(0i,02) plane, 03 = 0, and for a — 4, under which the original superpotential becomes 
harmonic. 

Recently, in the context of a three-field potential model was considered an hexagonal 
network of static classical configuration inside a topological soliton. Also, the 1/4 BPS 
equations of domain wall junction were first obtained by Gibbons and Townsend [21] and 
independently by Carrol et al. [22] . We point out that the superpotential model investigated 



23]. 



here can be applied to implement new string junctions by extended BPS domain walls 
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